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ABSTRACT 

The radii of young pre-main-sequence (PMS) stars in the Orion Nebula Cluster (ONC) 
have been estimated using their rotation periods and projected equatorial velocities. 
Stars at a given effective temperature have a spread in their geometrically estimated 
projected radii that is larger than can be accounted for with a coeval model, observa- 
tional uncertainties and randomly oriented rotation axes. It is shown that the required 
dispersion in radius (a factor of 2-3 full width half maximum) can be modelled in terms 
of a spread in stellar ages larger than the median age of the cluster, although the de- 
tailed star formation history cannot be uniquely determined using present data. This 
technique is relatively free from systematic uncertainties (binarity, extinction, variabil- 
ity, distance) that have hampered previous studies of the ONC star formation history 
using the conventional Hertzsprung-Russell diagram. However, the current ONC ro- 
tational data are biased against low luminosity objects, so the deduced dispersions 
in radius and inferred age are probably underestimates. In particular, the ages of a 
tail of PMS stars that appear to be ^ lOMyr old in the Hcrtzsprung-RusscU diagram 
cannot be verified with present data. If projected equatorial velocities were measured 
for these objects it could easily be checked whether their radii are correspondingly 
smaller than the bulk of the ONC population. 



Key words: 

M42 



stars: formation - methods: statistical - open clusters and associations: 



1 INTRODUCTION 

Whether star formation takes place rapidly on dynamical 
timescales or is a quasi-static process in which protostellar 
cores take many free-fall timescales to contract is one of the 
important debates in theoretical star formation (e.g. Shu, 
Adams & Lizano 1987; Hartmann, Ballesteros-Paredes & 
Bergin 2001; Tassis & Mouschovias 2004; Vazquez-Semadeni 
et al. 2005; Tan, Krumholz & McKee 2006). 

A key observational constraint on star formation 
timescales is the possibility of significant age spreads in 
young clusters and associations. Palla & Stahler (1999, 2000) 
have used ages estimated from the Hertzsprung-Russell (H- 
R) diagrams of several clusters to argue that star formation 
takes place over more than 10 Myr and at an accelerating 
rate towards the present day. Huff & Stahler (2006) show 
that this formation history appears to be similar in the in- 
ner and outer parts of the Orion Nebula cluster (ONC), 
suggesting that the acceleration was triggered by a global 
contraction of the parent cloud. However, these conclusions 
are not universally accepted. Hartmann (2001) has shown 
that errors in determining the luminosities of young stars, 
due to binarism, uncertain extinction, variability and accre- 



tion can lead to significant overestimates of the age spread 
in an H-R diagram and obscure the detailed star formation 
history. 

In this paper I explore an alternative technique for esti- 
mating the extent of luminosity and (inferred) age spreads in 
star forming regions using pre-main-sequence (PMS) stars in 
the ONC as an example. Rotational periods and projected 
equatorial velocities can be combined to give a geometric 
estimate of the stellar radius multiplied by the sine of an 
unknown rotation axis inclination, _Rsini. Assuming ran- 
dom axial orientations and using PMS evolutionary models, 
the distributions of Rsini can be modelled in terms of a 
dispersion in true radius and inferred age. This idea was 
stimulated by the work of Rhode, Herbst & Mathieu (2001), 
who calculated mean R sin i as a function of position in the 
H-R diagram, finding 3-sigma evidence that stars with lower 
luminosities did indeed have smaller radii. Here I provide a 
more sophisticated analysis of a larger dataset. This work 
follows closely on from the related analysis presented in Jef- 
fries (2007), where a similar approach was used to provide 
a revised distance estimate for the ONC. 

The data base of rotational measurements used in this 
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Figure 1. A Hertzsprung- Russell diagram for the 95 stars of the 
ONC rotation sample assuming an ONC distance of 392 pc (see 
text). Solid symbols indicate objects classified as accreting T- 
Tauri stars on the basis of an I — K colour excess greater than 
0.3 mag. Open symbols indicate non-accreting objects using the 
same criterion. Small symbols indicate all proper-motion selected 
members of the ONC from Hillenbrand (1997). The lines on the 
plot are isochrones from the models of D'Antona & Mazzitelli 
(1997). From the top downwards these isochrones are at 0.1, 0.3, 
1.0, 3.0, 10, 30 and 100 Myr respectively. 

paper is described in section 2. Section 3 presents the anal- 
ysis methods used in this work. The results are presented in 
section 4 and discussed in section 5. 



2 THE OBSERVATIONAL DATABASE 

The ONC is among the best studied star forming regions 
and the premier cluster for investigating star formation and 
early stellar evolution. It is relatively nearby, very young 
(392 ± 34 pc, ~l-2Myr - Jeffries 2007; 389 PC - Sand- 
strom et al. 2007) and contains a large population of stars 
and brown dwarfs covering the entire (sub) stellar mass spec- 
trum (0.01 < M/Mq < 30 - see Hillenbrand 1997; Slesnick, 
Hillenbrand & Carpenter 2004). 

Jeffries (2007) assembled a database of 95 K- and M- 
type PMS stars in the ONC, with model-dependent masses 
of about 0.2-2 Mq and which have proper-motion or radial 
velocity membership credentials, spectral types and lumi- 
nosities (Hillenbrand 1997) , effective temperatures based on 
the work by Cohen & Kuhi (1979), indicators of whether 
objects were actively accreting (based on I — K excess, 
from Hillenbrand [1997]), rotation periods (Herbst et al. 
2002) and usini (measured by Rhode et al. [2001] or Sicilia- 
Aguilar et al. [2005]). The typical precision of the measure- 
ments are ±1 spectral subclass, ±150 K in effective temper- 
ature, < 1% uncertainty in the periods and ~ 10% uncer- 
tainties in 17 sin i. The smallest «sini values that are taken 
to be reliably measured are 13.6 km s~^ and 10.0 km s~^ for 
the Rhode et al. and Sicilia-Aguilar et al. measurements re- 
spectively. 

There are various filters applied to this data and a num- 
ber of possible selection effects. The most important of these 
with regard to possible radius or age spreads are: (i) a bias 
against the inclusion of classical T-Tauri stars (CTTS), be- 



cause their rotation periods can be difficult to distinguish 
from other non-periodic variability; (ii) a bias against slowly 
rotating stars, either because their periods are too long to 
reliably detect, they are not so magnetically active or their 

V sin i is too small to be reliably measured; and (iii) the sam- 
ple of stars which are bright enough to obtain spectroscopy 
is biased towards optically brighter and hence intrinsically 
more luminous (and probably younger) stars. 

If CTTS are predominantly younger than non-accreting 
weak-lined T-Tauri stars (WTTS) then the first effect will 
bias the sample towards older ages. The second effect is more 
complicated. As PMS stars contract they spin-up, but they 
may also lose angular momentum through winds or star-disc 
interactions as they get older. The longer period stars that 
remain in our sample appear to have larger radii on aver- 
age (from radii estimated by Hillenbrand 1997) and may be 
younger. On the other hand, at a given rotation period a 
smaller (older) star will have a lower vsini. Therefore these 
selection effects probably bias our sample against the in- 
clusion of some of the youngest and oldest stars. Finally, 
any bias against low-luminosity stars will potentially ex- 
clude older stars, especially in the lowest mass objects in 
our sample. The discussion of completeness in Hillenbrand 
(1997) suggests that any incompleteness in the stars with 
spectroscopy will arise from the requirements of obtaining 

V sin i using high resolution spectroscopy, which is restricted 
to brighter stars. Spectral types can be obtained using low 
resolution spectroscopy for much fainter objects. 

An illustration of these potential effects is provided by 
comparing the H-R diagrams for our "rotation sample" and 
for proper-motion selected ONC members with tempera- 
tures and luminosities from Hillenbrand (1997). The H-R 
diagram is shown in Fig. 1, along with isochrones from 
DAntona & MazziteUi (1997). The luminosities from Hil- 
lenbrand (1997) have been adjusted to correspond to the 
392 pc found by Jeffries (2007), rather than the 470 pc as- 
sumed by Hillenbrand (1997). 

Model-dependent PMS stellar ages are estimated from 
these isochrones. Histograms of age are shown in Fig. [5] for 
the rotation sample and the full ONC sample, using lin- 
ear and logarithmic age bins. The median ages are 0.62 Myr 
and 1.03 Myr for the rotation sample and the full sam- 
ple respectively. The linearly binned age distribution of 
the Hillenbrand (1997) sample in Fig. [2] shows the pseudo- 
exponential decline characterised as "accelerating star for- 
mation" by Palla & Stabler (1999). This translates into a 
peak in the number of stars per logarithmic age interval 
at around 1 Myr. Fig. [5] also demonstrates that the tail of 
"older" (> 3 Myr) stars appears to be absent from the ro- 
tation sample. A Kolmogorov-Smirnov (K-S) test suggests 
that there is a < 1 per cent probability that the two age dis- 
tributions are consistent. This selection effect arises chiefiy 
from the bias towards brighter targets. Some of the age bias 
might be thought to arise because older, smaller objects 
could have unresolveable v sin i (see RebuU, Wolff & Strom 
2004). However, the H-R diagram-based age distribution for 
objects which have unresolved vsini but otherwise match 
the sample selection criteria is in fact similar to that of the 
rotation sample. 

Stars in the rotation sample are classified as CTTS or 
WTTS on the basis of whether their I — K colour excess is 
greater or less than 0.3 mag (see Hillenbrand et al. 1998). A 
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Figure 2. Deduced age distributions from Fig. [T]for tire rotation sample and for proper-motion selected ONC members from Hillenbrand 
(1997). The left panel is binned linearly showing the "exponentially accelerating" star formation rate proposed by Palla & Stahler (1999). 
The right panel is binned logarithmically. It is clear from this plot that a tail of older stars that is present among the proper-motion 
selected ONC members is not present in the rotation sample. 
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Figure 3. Deduced age distributions from the H-R diagram for the rotation sample and for proper-motion selected ONC members, but 
this time using the Siess et al. (2000) models to assign ages. Left and right panels show linearly and logarithmically binned ages. 



K-S test on the age distributions of these two populations 
reveals no significant evidence for a difference. 

The age distributions in Fig. [2] are of course model de- 
pendent. I have repeated these exercises using the Siess, Du- 
four & Forestini (2000) models (the Z = 0.02 variety) and 
these are shown in Fig. [S] The Siess et al. models give older 
ages and the evidence for "accelerating" star formation up 
to the present day now disappears to be replaced by a peak 
in the star formation rate in the recent past. However, the 
conclusions regarding a bias in the rotation sample and the 
missing tail of older stars remain. The median ages using 
the Siess et al. (2000) models are 1.55 Myr and 2.51 Myr for 
the rotation sample and full samples respectively. 



3 COMPUTATIONAL METHODS 
3.1 Monte Carlo Models 

The background to the method used in this paper is de- 
scribed by Jeffries (2007). The projected equatorial radius 
of a PMS star can be estimated from 



(i?sini)obs = — fobs (usin j)obs (1) 

where i is the inclination of the rotation axis to the line of 
sight, Pobs is the observed period, («sin j)obs is the observed 
projected equatorial velocity and k a constant dependent on 
the system of units in use. For any group of stars a mean 
i?sini could be found and divided by an average projection 
factor to obtain an estimate of the true radius. This was 
the approach adopted by Rhode et al. (2001), but is subject 
to a number of problems. First, many objects are not in- 
cluded in the considered samples because their v sin i values 
are too small to be resolved or their periods could not be 
measured because sini was too small. Second, uncertainties 
in the measured P and v sin i values themselves mean that 
the underlying distribution of projection factors is modified 
from a pure sinusoidal form. The aim here is to model the 
distribution of (J?sini)obs using a Monte Carlo simulation, 
starting from an assumed distribution of Rtme, the assump- 
tion of random axial inclination and taking account of the 
selection effects and measurement errors. 

If we accept that Pobs is related to Ptruo, the true ro- 
tation period in the absence of measurement uncertainties, 
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according to a normal error distribution characterised by a 
fractional measurement uncertainty 5P, then 

Pobs = Ptruc(l+(5p(7), (2) 

where U is a. random number drawn from a Gaussian distri- 
bution with mean of zero and unit standard deviation. The 
precision of the periods is generally very good and we as- 
sume 5p = 0.01 in all cases. Similarly, there is a relationship 
between the true and observed projected equatorial velocity 

(wsini)obs = "true sini (1 + S^U) , (3) 

where 5v is the fractional observational uncertainty in 
(«sini)obs (taken from the papers where the measurement 
was presented). In this case Vtme and sini are split into sep- 
arate factors in order to properly model the observational 
lower limits to v sin i and any selection effect on the possible 
values of i (see below). 

Equations 1-3 can be combined to give 

{Rsmi)ohs = Rtruc{l+SpUi){l + SvU2) sini, (4) 

where sini can be generated assuming random axial orien- 
tation and Ui and U2 are different random numbers taken 
from a unit Gaussian distribution. Note that in this formula- 
tion (_Rsini)obs is independent of the assumed Pt rue or ?^truc 
values, but that this is not the case once selection effects are 
considered (see below). 

For a given observational sample, 10'* randomised val- 
ues of (_Rsini)obs are generated for each star. The initial 
Rtrue distribution is taken from model isochrones at the TcS 
corresponding to each star in the observed sample. As a 
further refinement -Rtrue can be drawn randomly from a dis- 
tribution, either in terms of a spread about the isochronal 
radius, or specified as a distribution of age (see section . 

There are some further complications which are dealt 
with in the model. First, because the isochrones have a sig- 
nificant slope in the H-R diagram, then an uncertainty in T^s 
leads to additional dispersion in the expected Rtme- This is 
accounted for with a ±150 K Gaussian perturbation in Tcff 
when -Rtruo is calculated from the evolutionary models. 

Second, the observational sample will be incomplete for 
objects which have a low inclination angle because starspot 
modulation may be difficult to observe in these stars (see Jef- 
fries 2007 for a discussion). This is dealt with by adopting 
a threshold inclination angle ith, below which randomised 
trials will be rejected from the model distribution. This free 
parameter principally affects the low end of the R sin i dis- 
tribution. It is assumed that ith is independent of rotation 
period. 

There is also the issue of limits imposed by spectral res- 
olution on the values of {v sin i)obs that can be recorded. This 
is dealt with by discarding randomised trials which have a 

V sin i lower than the threshold appropriate for the dataset 
in question. The minimum measurable value of (wsini)obs 
is well defined, but the rejection process requires that an 
initial Vtma distribution is also specified. Here I follow the 
example in Jeffries (2007) and hypothesise a simple Vtmc 
distribution, which after multiplying by the projection fac- 
tor and accounting for observational uncertainties, provides 
a good match (measured with a K-S test) to the observed 

V sin i distribution. The exact choice of Vtmc distribution has 
very little effect on the simulated iisini distributions (see 
also Jeffries 2007). 



3.2 Age distributions 

A number of possible radius or age distributions can be 
tested to see whether they provide a reasonable description 
of any particular sample. The scenarios investigated are: (1) 
a co-eval population; (2) a Gaussian distribution of logj^Q R 
around an isochronal locus; (3) a radius calculated according 
to a Gaussian distribution in log^Q age; (4) a radius calcu- 
lated from an age distribution which exponentially decays 
beyond an initial starting age; and (5) a radius calculated 
according to the age distributions inferred from the tradi- 
tional H-R diagram (see Figs. [5] and 

Because i?true is expected to vary with T^s at a given 
age then the distributions of -Rtruc should be specified as 
a function of Teg. This problem is finessed by normalising 
J^truo by its value at an age of 3 Myr and at the T^s of a 
given star (for the observations) or trial (for the models). 
This effectively collapses the two-dimensional distribution 
along the isochrones, reducing the observed distribution to 
one-dimensional form and will be termed the distribution of 
normalised projected radii, i? sin i/iisMyr . 

The results from scenarios 1 and 2 will be largely in- 
dependent of the choice of evolutionary model as they rely 
only on normalising the i?truc values with an appropriate 
isochrone of R vs Tag. However, the absolute ages and age 
spreads in scenarios 3-5 could be very dependent on choice 
of evolutionary model. With this in mind the evolutionary 
models of D'Antona & Mazzitelh (1997 - hereafter DAM97) 
and the solar metallicity models of Siess, Dufour & Forestini 
(2000, hereafter SOO) have both been tested. 

There are too many free parameters and too few data 
points in the binned distribution of (_Rsini)obs (see sec- 
tion |4| to attempt an inversion to a true radius distribution 
or to perform chi-squared fitting - although this may be 
possible in future with more numerous data points. Instead, 
K-S tests are used to determine whether a model is capable 
of providing a satisfactory description of the observed (cu- 
mulative) (_Rsini)obs/i?3Myr distribution and to show which 
models can be ruled out by the data. In each scenario the 
parameter which determines the mean i? sin i/i?3Myr value, 
i.e. the central isochronal age, has been adjusted to min- 
imise the K-S statistic and maximise consistency between 
data and model. 

The main goals are to answer the following questions: 



(i) Is the normalised R sin i distribution of the ONC ro- 
tation sample consistent with a coeval sample that has no 
dispersion in radius about an isochronal value? 

(ii) If not, then what is the spread in radius about an 
isochronal value that could best explain the observed nor- 
malised i?sini distribution? 

(iii) If the dispersion in radius is modelled as an age dis- 
tribution then what is the age spread implied by the data 
and is any particular form of the age spread (Gaussian, ex- 
ponential) preferred? 

(iv) Is the distribution of normalised i?sini consistent 
with the age distribution inferred from the positions of the 
same stars in the H-R diagram? 
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Figure 4. Normalised {Rs'mi)^^,^ versus T^g for the rotation sample. (i?sini)oijs was calculated using equation [T] and then divided by 
the isochronal value of R at 3Myr. The two panels show the results of these calculations using the evolutionary models of DAM97 (left) 
and SOO (right). CTTS and WTTS (classified according to their near-IR excesses) are given different symbols. In each panel the solid 
lines indicate isochrones of normalised radius at 0.1, 0.3, 1.0, 3.0, 10 and 30Myr. 



4 RESULTS 

4.1 The projected stellar radii 

Figure |4] shows the normalised projected radii for the ro- 
tation sample as a function of T^s- Both the DAM97 and 
SOO models have been used to normalise the (i?sini)obs val- 
ues calculated from equation [1] and the results are shown 
in the left and right hand panels respectively, compared 
to isochrones of normalised radius from the same mod- 
els. In neither case is there any evidence for a trend in 
(7?sini)obs/^?3Myr with Tes (a least squares fit gives a 
flat line within 1-sigma uncertainty in both instances), 
which supports the approach of using the distribution 
of (7?sini)obs/i?3Myr as a Tcff-independent indicator of 
the true spread in radii (or age). The mean values of 
(i?sini)obs/7?3Myr are 1.51 ± 0.08 and 1.24 ± 0.07 for the 
DAM97 and SOO models respectively. Hence the mean age 
of the ONC (as judged by the stellar radii) is significantly 
younger according to the DAM97 models than according 
to the SOO models, but in both cases is < 3 Myr because 
(i?sini)obs should be a lower limit to -Rtrue- 

4.2 A coeval population 

The first distribution of -Rtruc tested was that of a coeval 
population. The age was specified and then the cumulative 
distribution of (/?sini)obs/^?3Myr tested against the Monte- 
Carlo predictions. The assumed age was adjusted to min- 
imise the K-S statistic and select the model i?truc distribu- 
tion least different to the observational data. This process 
was repeated using both the SOO and DAM97 evolutionary 
models to calculate i?3Myr. The numerical results are listed 
in Table [T] and illustrated in Fig [5] (the coeval models are 
the rows with a dispersion in radius [see section 4.3], set to 
ffr = 0). Initially I used Gaussian Tofj uncertainties of 150 K 
and a threshold inclination ith ~ 30° . The assumed vtmc dis- 
tribution consists of 20 per cent of stars following a uniform 
distribution for 10 < vtmc < 120kms~^ with the remaining 
80 per cent following an exponential distribution with a de- 
cay constant of 12kms~^ between 10 < Utrue < 120kms~^. 



Adopting this, the observed v sin i distribution is well repro- 
duced by the model with a K-S rejection probability of < 10 
per cent (see Jeffries 2007). 

For both the DAM97 and SOO models the data are in- 
consistent with the PMS population being drawn from a 
single coeval isochrone, at confidence levels of 99.6 and 96.3 
per cent respectively. These must be considered lower-limits 
to the model-rejection confidence level because the mean age 
was adjusted as a free parameter to minimise the K-S statis- 
tic. The "best-fitting" ages are log]^Q(age/Myr) of 5.88 and 
6.25 for the DAM97 and SOO models respectively. These are 
similar to, but slightly older than, the median ages found 
from the H-R diagrams in section 2, but note that values 
determined from the projected radii are independent of the 
assumed distance to the ONC and unaffected by the binary 
status of any of the stars. 

Fig. [5] shows that the coeval models are too narrow in 
Rsini. Obviously a radius (or age) spread could broaden 
the model distribution, but so too could a decrease in the 
assumed value of ith or an increase in the assumed TbH un- 
certainties. To test this I ran the coeval SOO models assum- 
ing ith ~ 15°. Table [T] shows the "best-fit" parameters and 
that the coeval model is still rejected with ~ 90 per cent 
confidence. Conversely, if ith were greater than 30° a coeval 
model would be an even poorer match to the data. I next 
tried broadening the model R sin i distribution by increasing 
the uncertainty in Teff. The model can still be rejected at 
^ 90 per cent confidence unless the Tch uncertainties are 
approximately doubled to 300 K, which is equivalent to a 
±2 spectral subclass (1-sigma) uncertainty in the spectral 
types. Note again that these confidence levels are lower lim- 
its because the mean age was tuned to minimise the K-S 
statistic. 

I conclude that unless the uncertainties in spectral clas- 
sification and Teff have been significantly underestimated by 
Hillenbrand (1997), then the observed values of iisini are 
inconsistent with coevality at a reasonably high level of con- 
fidence. 



© 2007 RAS, MNRAS OOO.ITlfTTI 



6 R.D. Jeffries 




SOO models 
log (age/Myr)=6.25 - 



2 3 
R sin(i)/R3 




DAM97 models 



log (age/Myr)=5.i 



2 3 4 

R sin(i)/R3 



E 

o 




2 3 
R sin(i)/R3 j^^^. 




R sin(i)/R3 



Figure 5. (Left) Normalised distributions of (_Rsini)obs/^3Myr (binned data) compared with models which assume a coeval population. 
The values of RsMyi taken from the models of Siess et al. (2000, top) or D'Antona & Mazzitelli (1997, bottom). (Right) A comparison 
of the cumulative distributions which are used for the K-S tests described in the text. These model distributions are too narrow to 
adequately represent the data, implying a spread in true radius. 



Table 1. A summary of the results obtained from comparing models with a dispersion in normalised radius (see text) with the observed 
distribution of i? sin j/i?3Myr ■ The comparisons were done using either the SOO or DAM97 models to calculate -RsMyr- Columns 1-4 
list the model parameters: the Gaussian sigma in logj^g(/?true/J?3Myr)i the threshold inclination for detection of periodic variability, the 
assumed T^g uncertainty and the mean logj^g age that best fits the data. Column 5 lists the K-S statistic -Dniax (see Press et al. 1992), 
column 6 gives the integrated probability that D > fmax, corresponding to one minus the model rejection confidence level, and column 
7 gives a comment on how well the model distribution matches the data. 



Ur ith ATcff Mean logig i'max P{D > Dmax) comment 
(deg) (K) Age/Myr 



Siess et al. (2000) models 



0.00 
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6.25 





144 


0.037 


too narrow 


0.00 


15 


150 


6.22 





125 


0.101 


marginally too narrow 


0.00 


30 


300 


6.25 





101 


0.275 


reasonable 


0.05 


30 


150 


6.26 





124 


0.103 


marginally too narrow 


0.15 


30 


150 


6.25 





041 


0.997 


good fit 


0.27 


30 


150 


6.26 





124 


0.103 


marginally too broad 


D'Antona & 


Mazzitelli (1997) models 










0.00 


30 


150 


5.88 





180 


0.004 


too narrow 


0.08 


30 


150 


5.89 





122 


0.111 


marginally to narrow 


0.15 


30 


150 


5.88 





060 


0.876 


good fit 


0.24 


30 


150 


5.85 





123 


0.108 


marginally too broad 
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Figure 6. The distribution of normalised (ijsin for ttie rotation sample compared with models that have a Gaussian distribution in 
logj^Q R around an isochronal value. The left hand panel shows the comparison for Gaussian sigmas of 0.05, 0.15 and 0.27 respectively - 
corresponding to the best fit value and the values below or above which the model can be rejected with > 90 per cent confidence. iJsMyr 
was calculated using the SOO models. The right hand panel shows the normalised probability distribution, integrated over all possible 
central Rtrue/RsMyr values (using either the SOO models [solid line] or the DAM97 models [dashed line]), that these distributions are 
consistent with the observed data as a function of the Gaussian dispersion, ar- The results using both evolutionary models suggest a 
significant and similar spread in radii. The vertical dashed lines enclose 90 per cent of the probability in the case of the SOO models. 



4.3 A radius dispersion 

Using model 2 (a Gaussian distribution of log^g R about a 
central value defined by a single isochronal age with a stan- 
dard deviation of cr^), I simulated a dispersion in radius 
about an isochronal value. The central isochronal age was 
again adjusted to minimise the K-S statistic and a range of 
ar was tested. The results, using the SOO models to calcu- 
late f?3Myr, are listed in Tableland illustrated in Fig. |6l 
If the dispersion increases to a,. > 0.05 then this model 
can no longer be rejected with > 90 per cent confidence. 
The "best fit" dispersion corresponds to ar = 0.15, whilst 
the model distribution becomes too broad and can be re- 
jected with > 90 per cent confidence if a,- > 0.27. The re- 
sults using the DAM97 models are similar - models with 
0.08 < ar < 0.24 cannot be rejected with better than 90 per 
cent confidence - but the required central isochronal age 
is still smaller (logjQ(age/Myr) is 5.88 compared with 6.25 
when using the SOO models - see Table [l|. 

The results above cannot be used to estimate confidence 
intervals on the best fitting value of ar, because the central 
isochronal age has been optimised in each case. The right 
hand panel of Fig. [6] shows a normalised probability distri- 
bution as a function of ar (calculated using the SOO models) , 
where the probability has been integrated over a wide range 
of possible mean ages. This illustrates that the best fitting 
ar is indeed about 0.15, with 90 per cent of the probabil- 
ity distribution contained between 0.08 and 0.24. The results 
using the DAM97 models are very similar and give a slightly 
narrower 90 per cent confidence interval of 0.10 < ar < 0.22. 

The Monte Carlo models suggest that the observed 
Rsini distribution can be well represented if the observed 
population has a spread, amounting to a full-width at half 
maximum of 0.2-0.5 dex, in the stellar radii around an 
isochronal value. This result is almost independent of which 
evolutionary models are considered. An obvious interpreta- 
tion is that this is caused by a spread in ages and this is 
investigated in the next subsection. 



4.4 An age spread 

Two separate models for an age spread were tested. The first 
was a Gaussian distribution of logj^g age around a central 
isochrone. The free parameters were the central age and the 
age dispersion (in logarithmic units), aa- The second model 
was an age distribution which is zero up to some starting 
age, jumps to a maximum and then decays exponentially 
with a decay constant, Aa, expressed in Myr. This latter 
model, with a suitably small starting age, represents the ex- 
ponentially accelerating star formation model advocated by 
Palla & Stahler (1999). The lowest starting age considered 
was 0.03 Myr. 

Ages drawn randomly from the Gaussian age distribu- 
tion were transformed into radii using the appropriate stellar 
models (at the T^s of each star in the observational dataset) 
and these radii were perturbed according to the 150 K Tcff 
uncertainties and then subjected to the measurement un- 
certainties, random axial orientations and selection effects 
before comparing the observed and modelled distribution of 
i?sini/i?3Myr. A grid of models covering a wide range of 
central ages and age dispersion was calculated for both the 
DAM97 and SOO models. Assuming that the "correct" solu- 
tion lay within this grid, I normalised the K-S probabilities 
and this gave a pair of relative probability grids which are 
shown in Figs. [7^ and[7lD. Projection of these grids onto the 
age dispersion axis gave a best fitting aa and a confidence 
interval (analogous to the procedure used to create Fig. [SJ 
which are listed in Table [2l The same procedure was used 
for the exponential age distribution, with the grid projected 
onto the Aa axis. The normalised probability grids for this 
distribution are shown in Figs.|8^ and|8)D and the numerical 
results are also listed in Tabled 

There are acceptable fits to the observational data for 
both Gaussian and exponential age distributions and for 
both the DAM97 and SOO models. Hence the current data 
are incapable of distinguishing between these possibilities. 
As expected, the 90 per cent confidence intervals on aa and 
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Figure 7. Images representing the relative probability of a Gaussian distribution of logjg age being able to describe the observed 
distribution of Rsini. The x-axes are the central isochronal age and the y-axes are the Gaussian dispersion (aa in dex). The contours 
enclose 68, 90 and 99 per cent of the probability distribution. The left hand plot corresponds to radii calculated using the DAM97 models; 
the right hand plot has radii calculated using the SOO models. Both plots have assumed ith = 30° (see Table [J} and both suggest a 
spread in age that is larger than the median age of the cluster. 



Table 2. A summary of the results of modelling the observed distribution of Rsini/ R^Myr with two assumed age distributions (see 
text). For each assumed distribution results are given for the best- fitting cases where radii wore calculated using either the DAM97 or 
SOO evolutionary models and for two values of the threshold inclination angle. The final column gives the 90 percent confidence interval 
for (Ta or Aa respectively. 



Gaussian dispersion in log age 



Model 


«th 


Log central a 


ge/Myr 


0"a 


^max 


P(D > Dmax) 


Ao-a 




(deg) 






(dex) 






(dox) 


DAM97 


30 


5.85 




0.48 


0.059 


0.886 


0.31-0.82 


DAM97 


15 


5.83 




0.43 


0.060 


0.873 


0.27-0.78 


SOO 


30 


6.22 




0.36 


0.045 


0.989 


0.20-0.56 


SOO 


15 


6.20 




0.34 


0.043 


0.994 


0.14-0.54 


Exponential age 


distribution 












Model 


«th 


Log zero point 


age/ Myr 




^max 


P(D > Dn,ax) 


AAa 




(deg) 






(Myr) 






(Myr) 


DAM97 


30 


4.55 




1.10 


0.058 


0.898 


0.72-1.30 


DAM97 


15 


4.95 




0.95 


0.057 


0.911 


0.58-1.15 


SOO 


30 


5.65 




1.90 


0.043 


0.993 


1.22-2.82 


SOO 


15 


5.65 




1.70 


0.037 


0.999 


0.98-2.59 



Aa do not encompass zero, which would correspond to a co- 
eval population. The dispersion in logj^Q age for the Gaussian 
distribution is quite similar for the DAM97 and SOO mod- 
els: both models suggest a full-width half maximum spread 
of factors of a few to > 10 and therefore an overall spread 
in absolute ages that is larger than the central age of the 
distribution and the median age of the ONC found from the 
H-R diagram. The exponential distribution has a favoured 



starting age of essentially zero for the DAM97 models and a 
decay time scale, Aa of 1.1 Myr. Again this indicates a spread 
of ages that is larger than the median age of the sample from 
the H-R diagram. The SOO models favour a non-zero maxi- 
mum in the age distribution, which is also seen in the H-R 
diagram-based age distribution shown in Fig. |3] although 
a very young starting age cannot be ruled out. The decay 
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Figure 8. Similar to Fig. [7]but for an age distribution that decays exponentially with time scale Xa (on the y-axis) from a maximum at 
a zero point age (on the x-axis) and is zero prior to this. Left and right hand plots correspond to radii calculated using the DAM97 and 
SOO models. Contours enclose 68, 90 and 99 per cent of the relative probability distributions. In both case Aa is similar to the median 
cluster age using the same evolutionary models. 



Table 3. A summary of the results of modelling the observed dis- 
tribution of /? sin i//?3Myr using the age distribution determined 
from the Hertzspring-Russell diagram. The column headings arc 
labelled as for Table 1. 



Model 


kh 
(dcg) 


^max 




comment 


DAM97 


30 


0.119 


0.130 


marginally narrow 


DAM97 


15 


0.089 


0.430 


good fit 


SOO 


30 


0.052 


0.957 


good fit 


SOO 


15 


0.066 


0.791 


good fit 



timescale is about 2 Myr and again this implies a spread of 
ages that is larger than the median age of the sample. 

I ran some more simulations using a smaller value of 
ith = 15°. The results for these are also given in Tabled 
As expected a smaller ith slightly reduces the age disper- 
sion required to explain the data, but does not affect the 
broad conclusions. Conversely, a larger ith would increase 
the required age dispersion. 



4.5 Ages based on the Hertzsprung-Russell 
diagram 

As a final test I constructed age distributions based upon 
the ages determined from the H-R diagram for the rota- 
tion sample, which are shown in the right hand panels of 
Figs.[2]and[3l These age distributions were used to construct 
model i?sini/i?3Myr distributions which in turn were com- 



pared with the observed i?sini/_R3Myr distributions via K-S 
tests. Initially, ith was fixed at 30°. A further set of models 
was generated with ith = 15° to investigate the influence of 
this parameter. 

The comparison between the model and observed distri- 
butions of i? sin i/_R3Myr are shown in Fig. [9] Numerical re- 
sults are listed in Table |3] All of the model distributions are 
formally compatible with the data and cannot be rejected 
at the 90 per cent confidence level, although the DAM97 
ith ~ 30° model barely passes this test. Hence the age distri- 
butions determined from the Hertzsprung-Russell diagram 
cannot be discounted as a valid description of the age dis- 
tribution required by the projected radii. Given the uncer- 
tainties in determining ages from the H-R diagram - bina- 
rity, accretion, reddening, variability (see Hartmann 2001) - 
that do not afflict determinations of the projected radii, this 
might seem somewhat surprising. One might have thought 
that the age distributions from the H-R diagram would be 
too broad because of these additional sources of error. 

The overall uncertainties in the luminosities have been 
considered by a number of authors (e.g. Hillenbrand 1997; 
Hartmann 2001; RebuU et al. 2004), who have concluded 
that Gaussian uncertainties in log^Q luminosity are in the 
range 0.16-0.20 dex. This leads approximately to Gaussian 
uncertainties in deduced age of 0.24-0.30 dex for PMS stars 
(Hartmann 2001). In section 4.4 I found that the age dis- 
persion required to explain the normalised projected radius 
distribution was probably larger than this. Hence it seems 
that the uncertainties that affect the H-R diagram-based 
ages are not large enough to significantly broaden the de- 
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Figure 9. Normalised projected radii distributions generated using the age distributions derived for the H-R diagrams in Fig. |2]and|3] 
The left and right hand panels show distributions using the DAM97 and SOO models respectively and that these match the observed 
distributions reasonably well. In each case models with ith = 15° and 30° are shown, demonstrating that this parameter does not affect 
the conclusions greatly. 



duced age distribution beyond that which is indicated by 
the more robust distribution of projected radii. 



5 DISCUSSION 

Palla & Stabler (1999) and Huff & Stabler (2006) presented 
conventional analyses of the ONC H-R diagram, conclud- 
ing that star formation began at a low level about 10 Myr 
ago and accelerated up to the present day. Slesnick et al. 
(2004) also found evidence for an older population in the 
ONC H-R diagram they assembled, but argued it might be 
due to younger stars which are viewed in scattered hght 
from circumstellar material. Hartmann (2001) asserted that 
these large age spreads could be an artefact of observational 
errors, binarity, variability, differential reddening and accre- 
tion luminosity that conspire to cause an observed spread in 
the H-R diagram that is not due to a significant spread in 
age. Hartmann (2001) demonstrated that symmetric uncer- 
tainties in logj^Q luminosity can lead to the appearance of an 
exponentially decaying linear age distribution (see Fig. [2]). 
Palla et al. (2007) support their argument for a true age 
spread by showing that a few percent of the ONC PMS 
population show evidence of lithium depletion which is con- 
sistent with ages > 10 Myr. 

The evidence presented in this paper supports the idea 
that at a given T^g, the ONC stars in our sample do ex- 
hibit a spread in radius (and hence luminosity) that is not 
consistent with isochronal radii at a single age. One inter- 
pretation of this is that the ONC stars have a range of ages. 
Our modelling shows that plausible age distributions can 
be hypothesised which do explain the observed dispersion 
in radii. These distributions have spreads in age which are 
larger than the median age of the sample. 

It is worth emphasising that the technique employed 
here is independent of whether stars have correctly esti- 
mated reddening, whether they are accreting, whether their 
distance is correctly estimated and other sources of un- 
certainty in calculating their luminosities (see Hillenbrand 
1997). The method does depend on reasonable estimates of 
the measurement uncertainties in vsini and Ted and the 



assumption that axial orientations are random. The latter 
assumption is difficult to support or falsify at present (see 
discussion in Jeffries 2007), although any concentration of 
inclination angle would require a greater spread of intrinsic 
radii or ages to match the observed spread in Rsini. 

A spread in ages is not the only way of explaining a 
dispersion in radius at a given TcS- For instance Tout, Livio 
& Bonnell (1999) show that an "age spread" in the H-R 
diagram could be produced in a coeval population of PMS 
stars if they have very different accretion histories. Irrespec- 
tive of whether the spread in radii found here is interpreted 
in terms of a spread in age, there are important consequences 
for those using the ONC (and other young star forming 
regions) to investigate evolutionary trends of angular mo- 
mentum loss, magnetic activity and circumstellar discs (e.g. 
Herbst et al. 2005) and for those using evolutionary models 
to estimate masses and find the (sub)stellar mass function. 
In these cases, the often-used assumption that stars in a 
given cluster are close-to-coeval or have similar properties 
(luminosity, radius) at a given mass or T^s, is not valid. 

The sample examined in this paper is subject to a num- 
ber of limitations. Foremost among these is that I have only 
been able to comment on the radius distribution of those 
stars for which periods and v sin i measurements are avail- 
able. Following the discussion in section 2, it is clear that 
this "rotation sample" is deficient in the small tail of stars 
which form the controversial elderly population discussed by 
Palla & Stabler (1999) and Palla et al. (2007). 

Whilst the H-R diagram-based age distribution does a 
reasonable job of reproducing the observed distribution of 
projected radii, it would be unwise to extrapolate and con- 
clude that the H-R diagram-based ages of the apparently 
elderly population are reliable. The likely non-Gaussian be- 
haviour of some of the causes of uncertainty in PMS lumi- 
nosities (e.g. accretion disk scattering or variability) means 
that this small fraction of outliers could yet be objects 
whose average intrinsic luminosity has been severely under- 
estimated. Fortunately, the determination of periods is still 
quite feasible for the fainter ONC members. Large telescopes 
with multi-object, high resolution spectroscopic capability 
should then be capable of determining v sin i for some frac- 
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tion of these. Even a small number of iisini values for these 
objects would be capable of revealing whether they have 
smaller radii than the rest of the ONC population, as ex- 
pected from their relative positions in the H-R diagram. 



6 SUMMARY 

I have taken a sample of PMS objects in the ONC with 
known rotation periods and projected equatorial velocities, 
and modelled the resultant distribution of projected equa- 
torial radii with a simple Monte-Carlo simulation. The sim- 
ulation assumes that rotation axes are randomly oriented in 
space and takes into account selection effects and observa- 
tional uncertainties in the data. A comparison between mod- 
els constructed using a variety of hypotheses concerning the 
intrinsic distribution of stellar radii leads to the following 
conclusions. 

(i) The observed distribution of projected equatorial radii 
is inconsistent with the hypothesis that the ONC sample 
is coeval with a high level of confidence. This conclusion 
can only be weakened if uncertainties in spectral types and 
inferred effective temperatures are significantly larger than 
estimated in the original literature sources. 

(ii) A spread in radius (assuming a Gaussian distribution 
in logj^o of 0.2-0.5 dex (full width half maximum) about 
an isochronal locus is required to match the data. This is 
almost independent of the choice of evolutionary model. 

(iii) If the radius dispersion is modelled in terms of a 
distribution of ages, then either a Gaussian or exponen- 
tially decaying form can feasibly reproduce the data. The 
absolute values of the median age and age spread or decay 
timescales depend on which evolutionary models are used, 
but age spreads must be greater than the median age of the 
ONC. 

(iv) The age distributions deduced from a conventional 
H-R diagram adequately reproduce the observed distribu- 
tion of projected radii, independent of which evolutionary 
models are used. This suggests that uncertainties in deriv- 
ing intrinsic luminosities are sufficiently small for the sample 
considered, that they do not mask the underlying spread in 
radii and inferred age. 

(v) The current data are biased in the sense that low- 
luminosity ONC members, and hence objects that are older 
in the H-R diagram, are not represented in published cat- 
alogues of projected equatorial velocity. This means that 
the radius and inferred age dispersions found here may be 
lower limits. In particular, the possibility that some clus- 
ter members are ^ 10 Myr older than the bulk of the ONC 
population cannot be tested. However, given projected equa- 
torial velocities for some of these objects it should easily be 
possible to check whether these objects have smaller radii 
consistent with their luminosity-based older ages. 



Hartmann L., Ballesteros-Paredes J., Bergin E. A., 2001, 
ApJ, 562, 852 

Herbst W., Bailer- Jones C. A. L., Mundt R., Meisenheimer 

K., Wackermann R., 2002, A&A, 396, 513 
Herbst W., Mundt R., 2005, Astrophys. J. , 633, 967 
Hillenbrand L. A., 1997, AJ, 113, 1733 
Hillenbrand L. A., Strom S. E., Calvet N., MerriU K. M., 

Gatley I., Makidon R. B., Meyer M. R., Skrutskie M. P., 

1998, Astron. J. , 116, 1816 
Huff E. M., Stabler S. W., 2006, ApJ, 644, 355 
Jeffries R. D., 2007, MNRAS, 376, 1109 
Palla P., Randich S., Pavlenko Y. V., Placcomio E., 

Pallavicini R., 2007, ApJ, 659, L41 
Palla P., Stahler S. W., 1999, ApJ, 525, 772 
Palla P., Stahler S. W., 2000, ApJ, 540, 255 
Press W. H., Teukolsky S. A., Vetterling W. T., Flannery 

B. P., 1992, Numerical Recipes in C: The art of scientific 
computing, 2nd ed.. Cambridge University Press, Cam- 
bridge, UK 

RebuU L. M., Wolff S. C, Strom S. E., 2004, AJ, 127, 1029 
Rhode K. L., Herbst W., Mathieu R. D., 2001, AJ, 122, 
3258 

Sandstrom K. M., Peek J. E. G., Bower G. C, Bolatto 
A. D., Plambeck R. L., 2007, ArXiv e-print 0706.2361, 
706 

Shu F. H., Adams P. C, Lizano S., 1987, ARA&A, 25, 23 
Sicilia-Aguilar A., Hartmann L. W., Hernandez J., Bricefio 

C, Calvet N., 2005, AJ, 130, 188 

Siess L., Dufour E., Forestini M., 2000, A&A, 358, 593 
Slesnick C. L., Hillenbrand L. A., Carpenter J. M., 2004, 

ApJ, 610, 1045 
Tan J. C, Krumholz M. R., McKee C. P., 2006, ApJ, 641, 

L121 

Tassis K., Mouschovias T. C, 2004, ApJ, 616, 283 
Tout C. A., Livio M., BonneU I. A., 1999, MNRAS, 310, 
360 

Vazquez-Semadeni E., Kim J., Shadmehri M., Ballesteros- 
Paredes J., 2005, ApJ, 618, 344 

This paper has been typeset from a TpjX/ BTjtX file prepared 
by the author. 



REFERENCES 

Cohen M., Kuhi L. V., 1979, ApJS, 41, 743 
D'Antona P., MazziteUi I., 1997, Mem. Soc. Astr. It., 68, 
807 

Hartmann L., 2001, AJ, 121, 1030 



© 2007 RAS, MNRAS OOO-ITlfTTI 



